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Abstract 



A search is presented for physics beyond the standard model (BSM) in final states 
with a pair of opposite-sign isolated leptons accompanied by jets and missing trans- 
verse energy. The search uses LHC data recorded at a center-of-mass energy a/s = 
7 TeV with the CMS detector, corresponding to an integrated luminosity of approxi- 
mately 5fb _1 . Two complementary search strategies are employed. The first probes 
models with a specific dilepton production mechanism that leads to a characteristic 
kinematic edge in the dilepton mass distribution. The second strategy probes models 
of dilepton production with heavy colored objects that decay to final states includ- 
ing invisible particles, leading to very large hadronic activity and missing transverse 
energy. No evidence for an event yield in excess of the standard model expectations 
is found. Upper limits on the BSM contributions to the signal regions are deduced 
from the results, which are used to exclude a region of the parameter space of the 
constrained minimal supersymmetric extension of the standard model. Additional 
information related to detector efficiencies and response is provided to allow testing 
specific models of BSM physics not considered in this paper. 
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1 Introduction 

In this paper we describe a search for physics beyond the standard model (BSM) in events 
containing a pair of opposite-sign leptons, jets, and missing transverse energy (E™ 1SS ), in a 
sample of proton-proton collisions at a center-of-mass energy of 7 TeV. The data sample was 
collected with the Compact Muon Solenoid (CMS) detector |T| at the Large Hadron Collider 
(LHC) in 2011 and corresponds to an integrated luminosity of 4.98 fb _1 . This is an update and 
extension of a previous analysis performed with a data sample of 34 pb _1 collected in 2010 0. 

The BSM signature in this search is motivated by three general considerations. First, new par- 
ticles predicted by BSM physics scenarios are expected to be heavy in most cases, since they 
have so far eluded detection. Second, BSM physics signals may be produced with large cross 
section via the strong interaction, resulting in significant hadronic activity. Third, astrophysical 
evidence for dark matter suggests (SHU that the mass of weakly-interacting massive particles is 
of the order of the electroweak symmetry breaking scale. Such particles, if produced in proton- 
proton collisions, could escape detection and give rise to an apparent imbalance in the event 
transverse energy. The analysis therefore focuses on the region of high E™ 1SS . An example of a 
specific BSM scenario is provided by R-parity conserving supersymmetric (SUSY) models, in 
which the colored squarks and gluinos are pair-produced and subsequently undergo cascade 
decays, producing jets and leptons flZHEI- These cascade decays may terminate in the produc- 
tion of the lightest SUSY particle (LSP), often the lightest neutralino, which escapes detection 
and results in large E™ 1SS . This LSP is a candidate for a dark matter weakly-interacting massive 
particle. Another BSM scenario which may lead to similar signatures is the model of universal 
extra dimensions (UED) 

The results reported in this paper are part of a broad program of BSM searches in events with 
jets and E™ ss , classified by the number and type of leptons in the final state. Here we describe a 
search for events containing an opposite-sign isolated lepton pair in addition to jets and E™ 1SS . 
We reconstruct electrons and muons, which provide a clean signature with low background. In 
addition, we reconstruct t leptons in their hadronic decay modes to improve the sensitivity to 
models with enhanced coupling to third generation particles. Complementary CMS searches 
with different final states have already been reported, for example in Refs. |IT0l[TT| . Results from 
the ATLAS collaboration in this final state using approximately 1-2 fb _1 have been reported in 

Refs. num. 

The analysis strategy is as follows. In order to select dilepton events, we use a preselec- 
tion based on that of the CMS top quark pair (tt) cross section measurement in the dilepton 
channel [14]; the details of this preselection are presented in Section [3] Reasonable agree- 
ment is found between the observed yields in data and the predictions from standard model 
(SM) Monte Carlo (MC) simulation. Two complementary search strategies are pursued, which 
are optimized for different experimental signatures. The first strategy is a search for a kine- 
matic edge fll5| in the dilepton (ee, fip) mass distribution. This is a characteristic feature 
of SUSY models in which the same-flavor opposite-sign leptons are produced via the decay 
x\ - > & where is the next-to-lightest neutralino, x\ is the lightest neutralino, 

and £ is a slepton. The second strategy is a search for an excess of events with dileptons accom- 
panied by very large hadronic activity and E™ 1SS . We perform counting experiments in four 
signal regions with requirements on these quantities to suppress the tt background, and com- 
pare the observed yields with the predictions from a background estimation technique based 
on data control samples, as well as with SM and BSM MC expectations. These two search ap- 
proaches are complementary, since the dilepton mass edge search is sensitive to new physics 
models that have lower E™ 1SS and hadronic energy, while the counting experiments do not as- 
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sume a specific dilepton production mechanism and are also sensitive to BSM scenarios that 
produce lepton pairs with uncorrelated flavor. 

No specific BSM physics scenario, e.g. a particular SUSY model, has been used to optimize 
the search regions. In order to illustrate the sensitivity of the search, a simplified and practi- 
cal model of SUSY breaking, the constrained minimal supersymmetric extension of the stan- 
dard model (CMSSM) I[l6j [TZB is used. The CMSSM is described by five parameters: the uni- 
versal scalar and gaugino mass parameters (mo and tni/2, respectively), the universal trilin- 
ear soft SUSY breaking parameter Aq, the ratio of the vacuum expectation values of the two 
Higgs doublets (tan and the sign of the Higgs mixing parameter ]i. Throughout the pa- 
per, four CMSSM parameter sets, referred to as LM1, LM3, LM6, and LM13 |[18|, are used 
to illustrate possible CMSSM yields. The parameter values defining LM1 (LM3, LM6, LM13) 
are m = 60 (330, 85, 270) GeV, m 1/2 = 250 (240,400,218) GeV, tan ,6 = 10 (20,10,40), A = 
(0,0, —553) GeV; all four parameter sets have }i > 0. These four scenarios are beyond the 
exclusion reach of previous searches performed at the Tevatron and LEP, and are chosen here 
because they produce events containing opposite-sign leptons and may lead to a kinematic 
edge in the dilepton mass distribution. These four scenarios serve as common benchmarks to 
facilitate comparisons of sensitivity among different analyses. 

2 The CMS detector 

The central feature of the CMS detector is a superconducting solenoid, 13 m in length and 
6 m in diameter, which provides an axial magnetic field of 3.8 T. Within the field volume 
are several particle detection systems. Charged particle trajectories are measured by silicon 
pixel and silicon strip trackers covering \tj\ < 2.5 in pseudorapidity where t] = — ln[tan0/2] 
with 9 the polar angle of the particle trajectory with respect to the counterclockwise proton 
beam direction. A crystal electromagnetic calorimeter and a brass /scintillator hadron calor- 
imeter surround the tracking volume, providing energy measurements of electrons, photons 
and hadronic jets. Muons are identified and measured in gas-ionization detectors embedded 
in the steel return yoke outside the solenoid. The detector is nearly hermetic, allowing energy 
balance measurements in the plane transverse to the beam direction. The first level of the CMS 
trigger system, composed of custom hardware processors, uses information from the calorime- 
ters and muon detectors to select, in less than 1 fis, the most interesting events. The High 
Level Trigger processor farm further decreases the event rate from around 100 kHz to around 
300 Hz, before data storage. Event reconstruction is performed with the particle-flow (PF) al- 
gorithm [19J, which is used to form a mutually exclusive collection of reconstructed particles 
(muons, electrons, photons, charged and neutral hadrons) by combining tracks and calorimeter 
clusters. A more detailed description of the CMS detector can be found elsewhere £Q. 

3 Event Selection 

The following samples of simulated events are used to guide the design of the analysis. These 
events are generated with either pythia 6.4.22 EDI , MadGraph 4.4.12 [21j, or POWHEG (221 
MC event generators using the CTEQ 6.6 parton density functions 1123 1 . The tt, W + jets, and 
W (V = W, Z) samples are generated with MadGraph, with parton showering simulated 
by pythia using the Z2 tune |24") . The single-top samples are generated with POWHEG. The 
Drell-Yan (DY) sample is generated using a mixture of MadGraph (for events with dilepton 
invariant mass above 50 GeV) and pythia (for events with dilepton invariant mass in the range 
10-50 GeV), and includes decays to the tt final state. The signal events are simulated using 
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pythia. The detector response in these samples is then simulated with a Geant4 model |25| 
of the CMS detector. The MC events are reconstructed and analyzed with the same software as 
is used to process collision data. Due to the varying instantaneous LHC luminosity the mean 
number of interactions in a single beam crossing increased over the course of the data-taking 
period to a maximum of about 15. In the MC simulation, multiple proton-proton interactions 
are simulated by PYTHIA and superimposed on the hard collision, and the simulated samples 
are reweighted to describe the distribution of reconstructed primary vertices in data [26|. The 
simulated sample yields are normalized to an integrated luminosity of 4.98 fb _1 using next-to- 
leading order (NLO) cross sections. 

Events in data are selected with a set of ee, ep., fiji, ex, and fir double-lepton triggers. Since the 
online reconstruction of hadronic-T decays (Th) is difficult, Th triggers are intrinsically prone to 
high rates. Therefore, for the analysis with two only, we use specialized triggers that rely on 
significant hadronic activity Hj, quantified by the scalar sum of online jet transverse energies 
with pj > 40 GeV, and EijF ss as well as the presence of two t^. The efficiencies for events 
containing two leptons passing the analysis selection to pass at least one of these triggers are 
measured to be approximately 1.00l{j;{g, 0.95 ± 0.02, 0.90 ± 0.02, 0.80 ± 0.05, 0.80 ± 0.05 and 
0.90 ± 0.05 for ee, eji, er^, yix^ and ThTh triggers, respectively. In the following, the simulated 
sample yields for the light lepton channels are weighted by these trigger efficiencies. For the Th 
channels the trigger simulation is applied to the MC simulation and then a correction is applied 
based on the measured data and MC efficiencies for these triggers. 

Because leptons produced in the decays of low-mass particles, such as hadrons containing b 
and c quarks, are nearly always inside jets, they can be suppressed by requiring the leptons to 
be isolated in space from other particles that carry a substantial amount of transverse momen- 
tum. The details of the lepton isolation measurement are given in Ref. |fl4|. In brief, a cone is 
constructed of size AR = a/ (A77) 2 + (A(p) 2 = 0.3 around the lepton momentum direction. The 
lepton relative isolation is then quantified by summing the transverse energy (as measured in 
the calorimeters) and the transverse momentum (as measured in the silicon tracker) of all ob- 
jects within this cone, excluding the lepton, and dividing by the lepton transverse momentum. 
The resulting quantity is required to be less than 0.15, rejecting the large background arising 
from QCD production of jets. 

The Th decays are reconstructed with the PF algorithm and identified with the hadrons-plus- 
strips (HPS) algorithm, which considers candidates with one or three charged pions and up to 
two neutral pions ||27| . As part of the Th identification procedure, loose isolation is applied for 
the ^ final states. Isolated electrons and muons can be misidentified as Th candidates. For this 
reason Th candidates are required to fail electron selections and not to match a muon signature 
in the muon system. 

Events with two opposite-sign isolated leptons are selected. At least one of the leptons must 
have pj > 20 GeV, both must have pj > 10 GeV, and the electrons (muons) must have \rj \ < 2.5 
(| 77 1 < 2.4). Electrons in the range 1.44 < \t]\ < 1.57 are excluded. In events containing a Th can- 
didate, both leptons must satisfy pj > 20 GeV and \ t]\ < 2.1, where the acceptance requirement 
is tightened so that the Th decay products are contained in the tracking detector in a manner 
that is consistent with the requirements of the triggers used for these events. In events with 
more than one opposite-sign pair that satisfy the selection requirements, the two oppositely- 
signed leptons with highest pj are chosen. Events with an ee or pp pair with invariant mass 
of the dilepton system between 76 GeV and 106 GeV or below 12 GeV are removed, in order to 
suppress Z/ 7* -4 11 events, as well as low-mass dilepton resonances. Events containing two 
electrons, two muons, or an electron and a muon are referred to as the "light lepton channels," 
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4 Search for a Kinematic Edge 



Table 1: Summary of event preselection requirements applied in the light lepton channels, 
hadronic-T channels, and the dilepton mass edge search of Section [4] The leading (trailing) 
lepton is the one with highest (second highest) pj. The requirements on jet multiplicity, scalar 
sum of jet transverse energies (Hj), missing transverse energy (£™ 1SS ), and dilepton mass are 
also indicated. 



Requirement 


light leptons 


hadronic-T 


edge search 


leading lepton 


e or }i, pj > 20 GeV 


e, ]i, or Th, pj > 20 GeV 


e or pj > 20 GeV 


trailing lepton 


e or ]i, pj > 10 GeV 


e, ]i, or Th, pj > 20 GeV 


e or ]i, pj > 10 GeV 


jet multiplicity 


"jets > 2 


"jets > 2 


"jets > 2 


Hj 


Hj > 100 GeV 


Hj > 100 GeV 


Hj > 300 GeV 


£miss 


Elf 55 > 50 GeV 


E™ iss > 100 GeV 


Ef ss > 150 GeV 


dilepton mass 


veto 76 < m ee , m w < 106 GeV 







while events with at least one Th are referred to as "hadronic-T channels." 

The PF objects are clustered to form jets using the anti-fry clustering algorithm [28] with the 
distance parameter of 0.5. We apply pj- and (/-dependent corrections to account for residual 
effects of nonuniform detector response, and impose quality criteria to reject jets that are consis- 
tent with anomalous detector noise. We require the presence of at least two jets with transverse 
momentum of pj > 30 GeV and \tj\ < 3.0, separated by AR > 0.4 from leptons passing the 
analysis selection. For each event the scalar sum of transverse energies of selected jets Hj must 
exceed 100 GeV. The E2p iss is defined as the magnitude of the vector sum of the transverse mo- 
menta of all PF objects, and we require E™ iss > 50 GeV (E™ ss > 100 GeV) in the light lepton 
(hadronic-T) channels. 

The event preselection requirements are summarized in Table [lj The data yields and corre- 
sponding MC predictions after this event preselection are given in Table [2] (light leptons) and 
Table |3] (hadronic-T). For the light lepton channels, the normalization of the simulated yields 
has been scaled based on studies of Z — > tl in data and in MC simulation, to account for effects 
of lepton selection and trigger efficiency and to match the integrated luminosity. As expected, 
the MC simulation predicts that the sample passing the preselection is dominated by lepton 
pair final states from tt decays (dilepton tt). The data yield is in good agreement with the 
prediction, within the systematic uncertainties of the integrated luminosity (2.2%) and tt cross 
section determination (12%) l29ti3TI . The yields for the LM1, LM3, LM6, and LM13 benchmark 
scenarios are also quoted. 

4 Search for a Kinematic Edge 

We search for a kinematic edge (end-point) in the dilepton mass distribution for same-flavor 
(SF) light-lepton events, i.e., ee or lepton pairs. Such an edge is a characteristic feature of, for 
example, SUSY scenarios in which the opposite-sign leptons are produced via the decay ^ — > 
££ — > X*(£ + £~- The model of UED can lead to a similar signature with different intermediate 
particles. In case of a discovery such a technique offers one of the best possibilities for model- 
independent constraints of the SUSY mass parameters [TBI . 

In contrast, for the dominant background tf as well as other SM processes such as WW and 
DY — > tt, the two lepton flavors are uncorrelated, and the rates for SF and opposite-flavor 
(OF) ey. lepton pairs are therefore the same. Hence we can search for new physics in the SF final 
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Table 2: Data yields and MC predictions in the light lepton channels after preselection, using 
the quoted NLO production cross sections a. The tt —> £ + £~ contribution corresponds to dilep- 
ton tt with no W — )• t decays, tt —¥ £ ± t t /t + t _ refers to dilepton tt with at least one W — > x 
decay, and tt — > ^ ± + jets/hadrons includes all other tt decay modes. The quoted cross sections 
for these processes include the relevant branching fractions. The LM points are benchmark 
SUSY scenarios, which are defined in the text. The MC uncertainties include the statistical 
component, the uncertainty in the integrated luminosity, and the dominant uncertainty from 
the tt cross-section determination. The data yield is in good agreement with the MC prediction, 
but the latter is not used explicitly in the search. The difference between the ee + jiji versus eji 
yields is due to the rejection of ee and ]i]i events with an invariant mass consistent with that of 
the Z boson. 



Sample 


<r[pb] 


ee 


W 


eji 


Total 


tt ->■ g+f 


7 


1466 ± 179 


1872 ± 228 


4262 ± 520 


7600 ± 927 


tt -> ^-rF/T+T" 


9 


303 ± 37 


398 ± 49 


889 ± 108 


1589 ± 194 


tt -> ^ ± + jets/hadrons 


141 


50 ± 6.2 


15 ± 1.9 


90 ± 11 


155 ± 19 




16677 


193 ±11 


237 ± 13 


312 ±15 


741 ± 26 


WW 


43 


55 ± 1.7 


66 ± 1.9 


151 ± 3.8 


272 ± 6.5 


wz 


18 


13 ± 0.4 


15 ± 0.4 


25 ± 0.6 


53 ± 1.3 


zz 


5.9 


2.6 ±0.1 


3.3 ± 0.1 


3.3 ±0.1 


9.1 ± 0.3 


Single top 


102 


95 ± 3.1 


120 ± 3.7 


278 ± 7.3 


492 ± 12 


W + jets 


96648 


47 ±11 


9.8 ± 4.6 


59 ± 12 


117 ± 16 


Total MC 




2224 ± 224 


2735 ± 281 


6069 ± 643 


11029 ± 1137 


Data 




2333 


2873 


6184 


11390 


LM1 


6.8 


272 ± 8.3 


342 ± 9.7 


166 ± 5.7 


780 ± 20 


LM3 


4.9 


107 ± 3.7 


125 ±4.1 


181 ± 5.5 


413 ± 11 


LM6 


0.4 


20 ± 0.6 


23 ± 0.7 


26 ± 0.8 


69 ± 1.7 


LM13 


9.8 


138 ± 6.6 


157 ± 7.0 


334 ± 12 


629 ± 19 



Table 3: Data yields and MC predictions in hadronic-T channels after preselection, using the 
quoted NLO production cross sections a. Diboson backgrounds comprise WW, WZ and ZZ 
events. The sum of simulated events is also split into events with a generated (MC, gen- 
uine Th) and events with a misidentified Th (MC, misidentified 7~h); the two contributions are 
equally important. The channel with two Th decays is not presented because the trigger is not 
efficient in the preselection region, due to the large Hj requirement. The uncertainty indicated 
represents both statistical and systematic components. 



Sample 


cr[pb] 


er h 




Total 


DY^ 11 


16677 


51 ± 12 


47 ±11 


98 ±22 


tt 


157.5 


165 ± 47 


205 ± 58 


370 ± 105 


Diboson 


66.9 


11 ± 2.0 


10.8 ± 1.9 


22 ± 3.6 


Single top 


102 


7.2 ± 2.6 


8.1 ± 2.7 


15 ± 4.8 


£ MC, genuine 




146 ± 39 


167 ± 44 


313 ± 83 


£ MC, misidentified T h 




89 ±24 


103 ± 27 


191 ± 51 


Total MC 


235 ± 62 


271 ± 72 


505 ± 134 


Data 


215 


302 


517 


LM1 


6.8 


36 ± 6.7 


46 ± 6.8 


82 ± 9.8 


LM3 


4.9 


28 ± 6.0 


18 ± 4.6 


46 ± 7.6 


LM6 


0.4 


2.8 ± 1.1 


4.2 ± 1.3 


7.0 ± 1.7 


LM13 


9.8 


90 ± 11 


118 ± 12 


208 ± 16 
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state and model the backgrounds using events in the OF final state. Thus the tt background 
shape is extracted from events with OF lepton pairs, and a fit is performed to the dilepton mass 
distribution in events with SF lepton pairs. 

In order to be sensitive to BSM physics over the full dilepton mass spectrum, events with a 
dilepton invariant mass trigg consistent with that of the Z boson are not rejected. This increases 
the DY contribution, which is compensated by an increase in the E^ ass > 150 GeV requirement 
(see Table [TJ. We then proceed to search for a kinematic edge in the signal region defined 
as Hj > 300 GeV. The invariant mass distributions of SF and OF lepton pairs are in good 
agreement with each other (see Fig. [I]). A fit is performed to the dilepton mass distribution 
with three candidate signal shapes, over a range of values on the position of the kinematic 
edge. 

The flavor-uncorrelated background, as a function of the invariant mass is parameterized 
as: 



where a m 1.4 describes the rising edge and b ~ 0.002 dominates the long exponential tail 
on the right hand side of the background shape; these parameters are extracted from the fit to 
data. 

We parametrize the signal shape with an edge model for two subsequent two-body decays, 
according to: 



For a = 1 this function describes a triangle convoluted with a Gaussian, which accounts for 
detector resolution effects. The resolution parameters for electrons a ee and muons a w are con- 



modelled by a Breit-Wigner function with the mass and width parameters fixed at the Z boson 
mass and width, convoluted with a Gaussian function to account for the detector resolution. 

We perform a simultaneous, extended, unbinned maximum likelihood fit to the distribution of 
dilepton mass for events containing ee, fiji (signal, DY and background model), and e}i pairs 
(background model only). The value of the kinematic edge position m max is varied, and the fit 
is performed for each value of this parameter. The shape parameters of the flavor-uncorrelated 
background that are free in the fit are assumed to be common in all categories, and the yields 
of signal (ftg), DY (mdy) an d background (n^) in these three categories are constrained using 
the ratio of muon to electron selection efficiencies R ?/e = 1.11 ± 0.05. This quantity is evaluated 
using studies of DY events in data and in MC simulation. 

The fit is performed in the signal region Hj > 300 GeV and E™ 1SS > 150 GeV. The SF events 
overlaid with the signal plus background fit, and the flavor-uncorrelated shape overlaid with 
OF events, are shown in Fig.JT] The results of the fit are displayed for a value of the kinematic 
edge position m max = 280 GeV, where the largest excess is observed. The local significance 
is 2.1(7" including statistical and systematic uncertainties. However, a correction for the look- 
elsewhere effect [32 1 reduces the global significance to 0.7a. The extracted signal yield includ- 
ing statistical uncertainty (ng = lllg'7) at this point is consistent with the background-only 
hypothesis, and we derive a 95% confidence level upper limit of n§ < 23 events for this kine- 
matic edge position. No evidence for a kinematic edge feature is observed in the dilepton mass 
distribution. 




(1) 




(2) 




is 
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m eH (GeV) m„ (GeV) 

Figure 1: Distribution of events (black points) and the results of the maximum likelihood fit 
(blue curve) to the dilepton mass distribution for events containing efi lepton pairs (left), and 
ee and jifi lepton pairs (right) in the signal region Hj > 300 GeV and Eij? iss > 150 GeV, that 
suppresses DY contributions almost completely. The signal hypothesis for a value of the kine- 
matic edge position m max = 280 GeV, corresponding to the largest local excess, is displayed. 
The shaded band represents the shape uncertainty of the background model. 

5 Counting Experiments 

We next proceed to search for an excess of events containing lepton pairs accompanied by large 
E™ 1SS and Hj. To look for possible BSM contributions, we define four signal regions that reject 
all but ~0.1% of the dilepton tf events, by adding the following requirements: 

• high-E™ 88 signal region : E™ 88 > 275 GeV, H T > 300 GeV, 

• high-H T signal region : E™ 88 > 200 GeV, H T > 600 GeV, 

• tight signal region : E™ 88 > 275 GeV, H T > 600 GeV, 

• low-H T signal region : E™ 88 > 275 GeV, 125 < H T < 300 GeV. 

The signal regions are indicated in Fig. [2] These signal regions are tighter than the one used in 
Ref. since with the larger data sample the tighter signal regions allow us to explore phase 
space farther from the core of the SM distributions. The observed and estimated yields in 
the high-E™ ss , high-Hx, and tight signal regions are used in the CMSSM exclusion limit in 
Section [7[ The low-Hj region has limited sensitivity to CMSSM models that tend to produce 
low-px leptons, since the large E™ 188 and low Hj requirements lead to the requirement of large 
dilepton pj. However, the results of this region are included to extend the sensitivity to other 
models that produce high-pj leptons. 

5.1 Light lepton channels 

The dominant background in the signal regions is dilepton tt production. This background 
is estimated using a technique based on data control samples, henceforth referred to as the 
dilepton transverse momentum (pj(££)) method. This method is based on the fact [33] that 
in dilepton tf events the pj distributions of the charged leptons (electrons and muons) and 
neutrinos are related, since each lepton-neutrino pair is produced in the two-body decay of the 
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W boson. This relation depends on the polarization of the W bosons, which is well understood 
in top quark decays in the SM ||34ti36l , and can therefore be reliably accounted for. In dilepton 
tf events, the values of pj(££) and the transverse momentum of the dineutrino system (pj(yv)) 
are approximately uncorrelated on an event-by-event basis. We thus use the observed pj(££) 
distribution to model the pj(vv) distribution, which is identified with E™ ss . Thus, we predict 
the background in a signal region S defined by requirements on E™ 1SS and Hj using the yield 
in a region S' defined by replacing the E™ 1SS requirement by the same requirement on pj(££). 

To suppress the DY contamination to the region S', we increase the E™ 1SS requirement to E™ 1SS > 
75 GeV for SF events and subtract off the small residual DY contribution using the R ou t/in tech- 
nique [14] based on control samples in data. This technique derives, from the observed DY 
yield in the Z mass region, the expected yield in the complementary region using the ratio 
Rout/in extracted from MC simulation. Two corrections are applied to the resulting prediction, 
following the same procedure as in Ref. |2J. The first correction accounts for the fact that we 
apply minimum requirements to E™ 1SS in the preselection but there is no corresponding re- 
quirement on pt(££). Since the E2p iss and pj(££) are approximately uncorrelated in individual 
dilepton tt events, the application of the E™ 1SS requirement decreases the normalization of the 
Pt{££) spectrum without significantly altering the shape. Hence, we apply correction factors 
K, which are extracted from data as K = 1.6 ± 0.1, 1.6 ± 0.4, 1.6 ± 0.4, and 1.9 ± 0.1 for the 
high-E™ lss , high-Hj, tight, and low-Hj signal regions, respectively. The uncertainty in K is 
dominated by the statistical component. The second correction factor Kq accounts for the W 
polarization in tt events, as well as detector effects such as hadronic energy scale; this correc- 
tion is extracted from MC and is K c = 1.6 ± 0.5, 1.4 ± 0.2, 1.7 ± 0.4, and 1.0 ± 0.4 for the four 
respective regions. The uncertainty in Kc is dominated by MC sample statistics and by the 
7.5% uncertainty in the hadronic energy scale in this analysis. 

Backgrounds from DY are estimated from data with the R ou t/in technique, which leads to an 
estimated DY contribution consistent with zero. Backgrounds from processes with two vector 
bosons as well as electroweak single top quark production are negligible compared with those 
from dilepton tt decays. 

Backgrounds in which one or both leptons do not originate from electroweak decays (misiden- 
tified leptons) are assessed using the "tight-to-loose" (TL) ratio (Rtl) method of Ref. |(T4l| . A 
misidentified lepton is a lepton candidate originating from within a jet, such as a lepton from 
semileptonic b or c decays, a muon from a pion or kaon decay-in-flight, a pion misidentified as 
an electron, or an unidentified photon conversion. The results of the tight-to-loose ratio method 
confirm the MC expectation that the misidentified lepton contribution is small compared to the 
dominant backgrounds. Estimates of the contributions to the signal region from QCD multijet 
events, with two misidentified leptons, and in W + jets, with one misidentified lepton in addi- 
tion to the lepton from the decay of the W, are derived separately. The contributions are found 
to be less than ~10% of the total background in the signal regions, which is comparable to the 
contribution in the control regions used to estimate the background from the pi{££) method. 
We therefore assign an additional systematic uncertainty of 10% on the background prediction 
from the pj{££) method due to misidentified leptons. 

As a validation of the pi{££) method in a region that is dominated by background, the pj(££) 
method is also applied in a control region by restricting Hj to be in the range 125-300 GeV. 
Here the predicted background yield is 95 ± 16 (stat) ± 40 (syst) events with E^ iss > 200 GeV, 
including the systematic uncertainties in the correction factors K and Kq, and the observed 
yield is 59 events. 

The data are displayed in the plane of E™ 1SS vs. Hj in Fig. The predicted and observed E™ ss 
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Figure 2: Distributions of E™ ss vs. Hj for data in the light lepton channels (left) and hadronic-T 
channels (right). The signal regions are indicated as hatched regions. The solid gray region is 
excluded at the preselection level. 



distributions are displayed in Fig. [3] A summary of these results is presented in Table|4j The SF 
and OF observed yields in the signal regions are quoted separately since many SUSY models 
lead to enhanced production of SF lepton pairs. For all signal regions, the observed yield is 
consistent with the predictions from MC and from the background estimate based on data. No 
evidence for BSM contributions to the signal regions is observed in the light lepton channels. 

5.2 Hadronic-T channels 

In the hadronic-T channels the background has two components of similar importance, events 
with a genuine lepton pair from dilepton tt production and events from semi-leptonic tt and 
W + jets production with a misidentified T^. Backgrounds are estimated separately with tech- 
niques based on data control samples. Other very small contributions from DY and diboson 
production with genuine lepton pairs ("MC irreducible") are estimated from simulation. 

The background with genuine lepton pairs is predicted by extending the pj(££) method. To 
translate the background prediction in the ee, eji, and channels into a prediction for the eTh, 
jlTh, and T^Th channels, a third correction factor is used. This correction, K T = 0.10 ± 0.01 for all 
signal regions, is estimated from simulation and accounts for the different lepton acceptances 
(~0.75), branching fractions (~0.56), and efficiencies (~0.24) in hadronic-T channels. This pro- 
cedure predicts the yield of the dilepton tt background with genuine hadronic t decays. 

The background with a reconstructed Th originating from a misidentified jet or a secondary 
decay is determined using a tight-to-loose ratio for t^ candidates measured in a dijet dominated 
data sample, defined as Hj > 200 GeV and E™ 1SS < 20 GeV. Tight candidates are defined 
as those that pass the full Th selection criteria. For the definition of loose candidates, the HPS 
isolation criterion is replaced by a looser requirement. The loose isolation requirement removes 
any Hj dependence of the tight-to-loose ratio; thus the measurement can be extrapolated to the 
signal regions. 

To determine the number of expected events including jets misidentified as t^ candidates in the 
signal region, the identification requirements for one t^ are loosened. The obtained yields are 
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5 Counting Experiments 
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Figure 3: The observed E™ ss distributions (red points) and E™ ss distributions predicted by 
the Pt(££) method (blue points with shaded uncertainty bands) in data for the region 125 < 
Hj < 300 GeV (upper left), H T > 300 GeV (upper right), and H T > 600 GeV (bottom). The 
uncertainty bands on the predicted E^ 1SS distribution are statistical, and also include systematic 
uncertainties for points in the signal regions, to the right of the vertical dashed line. The ratio 
of data to predicted background is also included. The error bars include the full uncertainties 
on the data and predicted background. 
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Table 4: Summary of results in the light lepton channels. The total SM MC expected yields 
(MC prediction), observed same-flavor (SF), opposite-flavor (OF), and total yields in the signal 
regions are indicated, as well as the predicted yields from the Pt(££) estimate. The the expected 
contributions from three benchmark SUSY scenarios are also quoted. The first uncertainty on 
the Pt(££) method prediction is statistical and the second is systematic; the systematic uncer- 
tainty is discussed in the text. The non-SM yield upper limit (UL) is a 95% CL upper limit on 
the signal contribution. 





high E™ 1SS 


high Hj 


tight 


low Hj 


MC prediction 


30 ± 1.2 


31 ±0.9 


12 ±0.6 


4.2 ± 0.3 


SF yield 


15 


11 


6 


3 


OF yield 


15 


18 


5 


3 


Total yield 


30 


29 


11 


6 


pj(££) prediction 


21 ± 8.9 ± 8.0 


22 ±7.5 ±6.9 


11 ±5.8 ±3.8 


12 ±4.9 ±5.7 


Observed UL 


26 


23 


11 


6.5 


Expected UL 


21 


19 


11 


8.6 


LM1 


221 ±5.1 


170 ± 4.5 


106 ± 3.5 


6.2 ±0.9 


LM3 


79 ± 2.4 


83 ± 2.5 


44 ±1.8 


2.3 ± 0.4 


LM6 


35 ± 0.6 


33 ± 0.5 


26 ± 0.5 


0.6 ±0.1 


LM13 


133 ± 5.5 


113 ±5.2 


65 ± 3.9 


4.1 ±0.9 



multiplied by the probability Pjl that a misidentified Th candidate passes the tight Th selection: 



1-Rtl(Pt,v)' 



A summation over Pjl evaluated for all Th candidates that pass the loose selection but not the 
tight selection gives the final background prediction in each signal region. 

The method is validated in tt simulation, where the agreement between the predicted and true 
yields is within 15%. We correct for a 5% bias observed in the simulation, and assign a 15% 
systematic uncertainty on the background prediction from the tight-to-loose ratio based on the 
agreement between prediction and observation in simulation and additional control samples 
in data. 

The results in the four signal regions are summarized in Table [5] The low-Hx region includes 
only ei"h and }it^ channels, because the T^Th trigger is inefficient in this region. In the high-E™ ss 
region the ThTh trigger is not fully efficient and an efficiency correction of 3% is applied to MC 
simulation. Good agreement between predicted and observed yields is observed. No evidence 
for BSM physics is observed in the hadronic-T channels. 

The results of observed yields and predicted backgrounds in all signal regions for different 
lepton categories are summarized in Fig. [4] 

6 Acceptance and Efficiency Systematic Uncertainties 

The acceptance and efficiency, as well as the systematic uncertainties in these quantities, de- 
pend on the process. For some of the individual uncertainties, it is reasonable to quote values 
based on SM control samples with kinematic properties similar to the SUSY benchmark mod- 
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6 Acceptance and Efficiency Systematic Uncertainties 



Table 5: Summary of the observed and predicted yields in the four signal regions for hadronic- 
t channels. The first indicated error is statistical and the second is systematic; the systematic 
uncertainties on the Rtl ratio and pj(££) method predictions are discussed in the text. The 
non-SM yield upper limit is a 95% CL upper limit on the signal contribution in each signal 
region. 





Vno-Vi F mi SS 
lllgll iZ*y 


lllgll ITT 


tignt 


iuw nj 


MC, genuine 


5.8 ±2.3 


3.7 ±1.6 


2.0 ±1.2 


0.4 ± 0.2 


MC, misidentified Th 


1.4 ±0.5 


2.8 ±1.3 


0.2 ±0.1 


0.2 ±0.1 


Total MC 


7.1 ±2.5 


6.5 ±2.3 


2.2 ± 1.2 


0.7 ±0.3 


pj{t£) prediction 


2.1 ±0.9 ±0.8 


2.2 ±0.8 ±0.9 


1.1 ±0.6 ±0.4 


1.2 ±0.5 ±0.4 


Rtl prediction 


5.1 ±1.7 ±0.8 


3.6 ±1.4 ±0.5 


2.7 ±1.3 ±0.4 


< 0.9@95%CL 


MC irreducible 


1.3 ±0.5 ±0.2 


0.7 ±0.3 ±0.1 


0.2 ±0.1 ±0.1 


0.1 ±0.1 ±0.1 


yj predictions 


8.5 ±2.0 ±1.1 


6.5 ±1.6 ±1.0 


4.0 ±1.4 ±0.6 


1.3 ±0.5 ±0.5 


Total yield 


8 


5 


1 





Observed UL 


7.9 


6.2 


3.7 


3.1 


Expected UL 


8.1 


7.2 


5.7 


3.9 


LM1 


32 ±11 


14 ±6.1 


8.1 ±4.2 




LM3 


11 ±4.2 


11 ±5.1 


8.0 ±4.9 




LM6 


4.5 ±1.5 


5.1 ±1.6 


4.2 ±1.6 


0.4 ± 0.4 


LM13 


69 ±17 


52 ± 8.2 


39 ± 9.8 
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Figure 4: Summary of the background predictions from tight-to-loose ratio, px(^)-method and 
MC, and observed yields in the signal regions. 
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els. For others that depend strongly on the kinematic properties of the event, the systematic 
uncertainties must be quoted model-by-model. 

The systematic uncertainty in the lepton acceptance consists of two parts: the trigger efficiency 
uncertainty and the identification and isolation uncertainty. The trigger efficiency for two lep- 
tons of pi > 10 GeV, with one lepton of pj > 20 GeV is measured using samples of Z — > ££, 
with an uncertainty of 2%. The simulated events reproduce the lepton identification and isola- 
tion efficiencies measured in data using samples of Z — > 11 within 2% for lepton pj > 15 GeV 
and within 7% (5%) for electrons (muons) in the range pj = 10-15 GeV. The uncertainty of the 
trigger efficiency (5%) of the triggers is estimated with the tag-and-probe method [37|. The 
Th identification efficiency uncertainty is estimated to be 6% from an independent study using 
a tag-and-probe technique on Z — > tt events. This is further validated by obtaining a Z — > tt 
enhanced region showing consistency between simulation and data. Another significant source 
of systematic uncertainty is associated with the jet and E™ ss energy scale. The impact of this 
uncertainty is final-state dependent. Final states characterized by very large hadronic activity 
and E™ 1SS are less sensitive than final states where the E™ 1SS and Hj are typically close to the 
minimum requirements applied to these quantities. To be more quantitative, we have used the 
method of Ref. |T4| to evaluate the systematic uncertainties in the acceptance for three bench- 
mark SUSY points. The energies of jets in this analysis are known to within 7.5%; the correction 
accounting for the small difference between the hadronic energy scales in data and MC is not 
applied |38| . 

The uncertainty on the LM1 signal efficiency in the region Hj > 300 GeV, E™ 1SS > 150 GeV 
used to search for the kinematic edge is 6%. The uncertainties for the four benchmark SUSY 
scenarios in the signal regions used for the counting experiments of Section [5] are displayed in 
Table [6] The uncertainty in the integrated luminosity is 2.2%. 

Table 6: Summary of the relative uncertainties in the signal efficiency due to the jet and E™ ss 
scale, for the four benchmark SUSY scenarios in the signal regions used for the counting exper- 
iments of Section IBl 



Signal Model 


high E™ 88 


high Hj 


tight 


low Hj 


LM1 


22% 


33% 


40% 


19% 


LM3 


26% 


34% 


42% 


18% 


LM6 


11% 


15% 


19% 


10% 


LM13 


26% 


31% 


40% 


14% 



7 Limits on New Physics 
7.1 Search for a kinematic edge 

An upper limit on the signal yield is extracted from the fit to the dilepton mass distribution, 
assuming the triangular shape (a = 1) of Eq. ((2). The 95% CL upper limit is extracted us- 
ing a hybrid frequentist-bayesian CLg method [39], including uncertainties in the background 
model, resolution model and Z-boson yield. We scan the position of the kinematic edge m max 
and extract a signal yield upper limit for each value, as shown in Fig. [5] The extracted up- 
per limits on n§ vary in the range 5-30 events; these upper limits do not depend strongly on 
the choice of signal shape parameter when using two different shapes specified by a concave 
(a = 4) and convex curvature (hatched band). 
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7 Limits on New Physics 
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Figure 5: A CLg 95% CL upper limit on the signal yield ftg as a function of the endpoint in 
the invariant mass spectrum m max assuming a triangular shaped signal (black dots and thick 
line). The hatched band shows the variation of the expected limit (thin line) assuming two 
alternate signal shapes, with the alternative expected limits corresponding to the boundary of 
the hatched band. The SUSY benchmark scenarios LM1, LM3 and LM6 are shown with their 
expected yields and theoretical positions of the corresponding kinematic dilepton mass edges. 
The LM1 (LM3) yield is scaled to 20% (40%) of its nominal yield. At LM3 and LM6 a three-body 
decay is present; thus the shape of the kinematic edge is only approximately triangular. 

7.2 Search for an excess of events with large E!f iss and H T 

In this section we use the results of the search for events with light leptons accompanied by 
large E™ 1SS and Hj reported in Section [5] to exclude a region of the CMSSM parameter space. 
The exclusion is performed using multiple, exclusive signal regions based on the high-E^ lss , 
high-Hx, and tight signal regions, divided into three non-overlapping regions in the E™ ss vs. 
Hj plane. The results are further divided between the SF and OF final states in order to im- 
prove the sensitivity to models with correlated dilepton production leading to an excess of SF 
events, yielding a total of six signal bins, as summarized in Table |7| The use of multiple, dis- 
joint signal regions improves the sensitivity of this analysis to a specific BSM scenario. The 
predicted backgrounds in the SF and OF final states are both equal to half of the total predicted 
background, because the tt events produce equal SF and OF yields. The inputs to the upper 
limit calculation are the expected background yields and uncertainties from the pj(££) method, 
the expected signal yields and uncertainties from MC simulation, and the observed data yields 
in these six regions. The exclusion is performed with the CLg method. In the presence of a 
signal, the pj(££) background estimate increases due to signal events populating the control 
regions. To correct for this effect, for each point in the CMSSM parameter space this expected 
increase is subtracted from the signal yields in our search regions. 

The SUSY particle spectrum is calculated using SoftSUSY f40l| , and the signal events are gen- 
erated at leading order (LO) with pythia 6.4.22. We use NLO cross sections, obtained with the 
program PROSPINO Elf . Experimental uncertainties from luminosity, trigger efficiency, and 
lepton selection efficiency are constant across the CMSSM plane, while the uncertainty from 
the hadronic energy scale is assessed separately at each CMSSM point taking into account the 
bin-to-bin migration of signal events. The variation in the observed and expected limits due to 
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Table 7: Summary of results in the light lepton channels used for the CMSSM exclusion of 
Section [7] Details are the same as in Table [4] except that these results are divided into three 
non-overlapping regions defined by E^ iss > 275 GeV, H T 300-600 GeV (SRI), E^ iss > 275 GeV, 
H T > 600 GeV (SR2, same as the "tight" signal region), and E™ iss 200-275 GeV, Hj > 600 GeV 
(SR3). The regions are further divided between same-flavor (SF) and opposite-flavor (OF) lep- 
ton pairs. 





SRI 


SR2 


SR3 


SF yield 


9 


6 


5 


OF yield 


10 


5 


13 


pj(££) prediction 


5.7 ±5.1 ±2.8 


5.3 ±4.1 ±1.9 


5.6 ±3.4 ±2.1 



the theoretical uncertainties, including renormalization and factorization scale, parton density 
functions (PDFs), and the strong coupling strength as ff42"J , are indicated in Fig. [6] as separate 
exclusion contours. These results significantly extend the sensitivity of our previous results f2~|. 
The LEP-excluded regions are also indicated; these are based on searches for sleptons and 
charginos fl43| . 

8 Additional Information for Model Testing 

Other models of new physics in the dilepton final state can be constrained in an approximate 
way by simple generator-level studies that compare the expected number of events in the data 
sample corresponding to an integrated luminosity of 4.98 fb _1 with the upper limits from Sec- 
tion [7] The key ingredients of such studies are the kinematic requirements described in this 
paper, the lepton efficiencies, and the detector responses for Hj and E™ 1SS . The trigger efficien- 
cies for events containing ee, e/Y or pfi lepton pairs are 100%, 95%, and 90%, respectively. For 
eih, and pr^ the efficiency is ~80% 11371 . The trigger used for t^tj^ final states has an efficiency 
of 90%. 

We evaluate the light lepton, hadronic-T, E™ ss , and Hj selection efficiencies using the LM6 
benchmark model, but these efficiencies do not depend strongly on the choice of model. Jets at 
the generator-level are approximated as quarks or gluons produced prior to the parton show- 
ering step satisfying pj > 30 GeV and \t]\ < 3. Generator-level leptons are required to satisfy 
px > 10 GeV and \t]\ < 2.5 and not to overlap with a generator-level jet within AR < 0.4. For 
generator level Th the visible decay products are required to satisfy the tighter pj > 20 GeV 
and 1 1]\ < 2.1 selection. The generator-level E™ 1SS is the absolute value of the vector sum of the 
transverse momenta of invisible particles, e.g., neutrinos and lightest supersymmetric parti- 
cles. The lepton selection efficiencies as a function of generator-level pj are displayed in Fig. [Z| 
The efficiency dependence can be parameterized as a function of pj as 

f{p T ) =e o{erf[(p T -C)/ t 7]}±e c {l-erf[(p T -C)/ t 7]}, (3) 

where erf indicates the error function, £oo gives the value of the efficiency plateau at high mo- 
menta, C is equal to 10 GeV, ec gives the value of the efficiency at pj = C, and c describes how 
fast the transition is. The parameterization is summarized in Table |H]for electrons, muons, and 
taus. 

The E^ 1SS and Hj selection efficiencies are displayed in Fig.|8]as a function of the generator-level 
quantities. These efficiencies are parameterized using the function: 
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8 Additional Information for Model Testing 
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Figure 6: The observed 95% CL exclusion contour (solid thick red line), the expected exclusion 
contour (solid thin blue line), the variation in the observed exclusion from the variation of PDF, 
renormalization and factorization scales, and a$ theoretical uncertainties (dashed red lines), 
the ±1(7 uncertainty in the median expected exclusion (dotted blue lines), and the observed 
exclusion contour based on 34 pb _1 2010 data in the opposite-sign dilepton channel (dark blue 
shaded region), in the CMSSM (mo,m 1/2 ) plane for tan/3 = 10, Aq = OGeV and ]i > 0. The 
area below the red curve is excluded by this search. Exclusion limits obtained from the LEP ex- 
periments are presented as shaded areas in the plot. The thin grey lines correspond to constant 
squark and gluino masses. The LM benchmark SUSY scenarios are also indicated. The LM3 
and LM13 benchmark scenarios have values of tan /3 and/ or Aq that differ from 10 and GeV, 
respectively, but both are also excluded by the results of this search; see the text of Section[l]for 
the full definitions of these scenarios. 



Table 8: Values of the fitted parameters in Eq. (pi for the lepton selection efficiencies of Fig. 7 



Parameter 


e 






C 


10 GeV 


10 GeV 


10 GeV 




0.78 


0.89 


0.44 




0.34 


0.62 


0.31 


a 


18 GeV 


30 GeV 


13 GeV 
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Figure 7: The efficiency to pass the light lepton (left), and hadronic-r (right) selection as a 
function of the generator-level pj (visible pj). These efficiencies are calculated using the 
LM6 MC benchmark. 
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Figure 8: The efficiency to pass the signal region E™ 1SS (left), and Hj (right) requirements 
as a function of the generator-level quantities. The vertical lines represent the requirements 
applied to the reconstruction-level quantities. These efficiencies are calculated using the LM6 
MC benchmark, but they do not depend strongly on the underlying physics. 
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9 Summary 



/(*) = ^(erf((*-C)Ar) + l), (4) 

where 6oo gives the value of the efficiency plateau at high x, C is the value of x at which the 
efficiency is equal to 50%, and a describes how fast the transition is. The values of the fitted 
parameters are quoted in Table [9] 

Table 9: Values of the fitted parameters in Eq. jij) for the E™ ss and Hj selection efficiencies of 
Fig-0 



Parameter 


E™ 1SS > 150 GeV 


E™ ss > 200 GeV 


E™ 1SS > 275 GeV 




1.00 


1.00 


1.00 


C 


157 GeV 


211 GeV 


291 GeV 


a 


33 GeV 


37 GeV 


39 GeV 


Parameter 


H T > 125 GeV 


Hj > 300 GeV 


Hj > 600 GeV 




1.00 


1.00 


0.99 


c 


124 GeV 


283 GeV 


582 GeV 


a 


56 GeV 


75 GeV 


93 GeV 



This efficiency model has been validated by comparing the yields from the full reconstruction 
with the expected yields using generator-level information only and the efficiencies quoted 
above. In addition to the LM1, LM3, LM6 and LM13 benchmarks considered throughout this 
paper, we have tested several additional benchmarks (LM2, LM4, LM5, LM7, and LM8) |18|. 
In general we observe agreement between full reconstruction and the efficiency model within 
approximately 15%. 

9 Summary 

We have presented a search for physics beyond the standard model in the opposite-sign dilep- 
ton final state using a data sample of proton-proton collisions at a center-of-mass energy of 
7TeV. The data sample corresponds to an integrated luminosity of 4.98 fb _1 , and was collected 
with the CMS detector in 2011. Two complementary search strategies have been performed. 
The first focuses on models with a specific dilepton production mechanism leading to a char- 
acteristic kinematic edge in the dilepton mass distribution, and the second focuses on dilepton 
events accompanied by large missing transverse energy and significant hadronic activity. This 
work is motivated by many models of BSM physics, such as supersymmetric models or mod- 
els with universal extra dimensions. In the absence of evidence for BSM physics, we set upper 
limits on the BSM contributions to yields in the signal regions. Additional information has 
been provided to allow testing whether specific models of new physics are excluded by these 
results. The presented result is the most stringent limit to date from the opposite-sign dilepton 
final state accompanied by large missing transverse energy and hadronic activity. 
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